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Abstract 

We examine the effect of a prior that favours low values of fine-tuning on Bayesian 
multi-dimensional fits of the constrained minimal supersymmetric standard model 
(CMSSM or mSUGRA) to current data. The dark matter relic density, the anoma- 
lous magnetic moment of the muon and the branching ratio of b — > 57 are all used 
to constrain the model via a Markov Chain Monte Carlo sampler. As a result of 
the naturalness prior, posterior probability distributions skew towards lighter higgs 
and sparticle masses, the effect being most pronounced in the gaugino sector. In- 
terestingly, slepton masses are an exception and skew towards heavier masses. The 
lightest CP-even Higgs /i°-pole annihilation mechanism becomes allowed at the 2a 
level for the latest combination of measurements of mt = 172.7 ±2.9 GeV, provided 
we allow for a theoretical error in the prediction of its mass m^o . m h o is constrained 
to be less than 120 GeV at the 95% C.L. Probing the branching ratio of B s — > 
to the level of 2xl0 -8 , as might be achieved by the Tevatron experiments, would 
cover 32% of the probability density, irrespective of which of the two priors is used. 



1 Introduction 



Although weak-scale supersymmetry (SUSY) [1] arguably offers the best can- 
didate for physics beyond the Standard Model, it potentially suffers from a 
fine-tuning problem [2]. Radiative electroweak symmetry breaking conditions 
imply that the mass of the Z°-boson Mz is related to the superpotential [i- 
term, the two soft supersymmetry breaking Higgs mass parameters 2 an d 
the ratio of the Higgs vacuum expectation values V2/V1 = tan/3 by 

M z _ m 2 Hl - m 2 H2 tan 2 g 2 
2 _ tan 2 /3-l [i) 
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at tree level [1]. In the following numerical calculations, the full one-loop 
MSSM corrections are added to this relation. Experiments at LEP2 and the 
Tevatron have placed increasingly severe bounds upon sparticle and Higgs 
masses [2,3] and in constrained models they push up the minimum value of fj, 
to several hundred GeV. In order to reproduce the observed value of M z = 91 
GeV, there then must be a large degree of cancellation between the a priori 
unrelated terms on the right hand side of Eq. 1. This is considered to be un- 
natural. In order to quantify the amount of fine-tuning, one can construct [2] 
a fine-tuning measure q for a parameter of the model Pi by noting that M z 
is unstable to small changes of Pi as a result of the fine-tuning: 



d\nM z 
d In pi 



c = maxjcj}. 



(2) 



q then measures the fractional change in M z (the partial derivative is calcu- 
lated while not using the constraint in Eq. 1) produced by some small fractional 
change in the parameter p^. We shall here use the largest q to quantify the 
amount of fine-tuning c for a given point in parameter space. 

It is our purpose in this paper to study the effect of a prior probability dis- 
tribution that favours a low value of c in combined fits to the CMSSM. We 
wish to extract any collider implications, as well as implications for the dark 
matter annihilation mechanism. The MSSM contains several good candidates 
for cold dark matter, the most studied being the lightest neutralino Xv The 
assumption that thermally produced neutralinos make up the dominant com- 
ponent of dark matter leads to strict constraints upon the parameter space 
of the MSSM. We shall here take this assumption, as well as the assumption 
of R-parity (which provides stability for the lightest neutralinos). In order to 
keep the dimensionality of the parameter space down, we shall further spe- 
cialise to the CMSSM. The CMSSM unifies all SUSY breaking scalar mass 
parameters to m , the Standard Model gaugino mass parameters to My 2 and 
the SUSY breaking trilinear scalar couplings to A at a high energy scale, 
commonly taken to be the scale of electroweak gauge unification Mqut- Our 
choice of parameters from which to calculate the fine-tuning includes the SUSY 
breaking Higgs bilinear B: 

Pi E {m , A , Mi/2, B, /i}. (3) 



B is used in preference to tan/5 since it directly appears in the Lagrangian 
and is therefore considered to be more fundamental. Our definition of c does 
not include the sensitivity to the GUT-scale top Yukawa coupling h t . We can 
argue about whether to include it or not in the definition since its value is 
determined by aspects of flavour physics in the model, an orthogonal aspect 
to the SUSY breaking aspects studied in the present paper. Were we to include 
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it in the definition of c, a strong suppression of the high tuq regime would be 
the likely result of the naturalness prior, since Ch t is very high there [4]. 

In ref. [5], exclusion limits from LEP1 and LEP2 were placed upon the CMSSM 
and a probability of 1 /c was applied to the remaining parameter space. This 
resulted in probability distributions for sparticle masses and various observ- 
ables. This work is similar in spirit to the ingredient of the naturalness prior 
that we will follow. A more common approach in the literature involves ap- 
plying separate 95% C.L. constraints upon a two-dimensional sub-space of 
the model, see recent examples in ref. [6]. The advantage of this simple ap- 
proach is that it is easy to discern which constraints act on different parts of 
parameter space. On the other hand, it was argued in Ref. [7] that a more 
complete approach would perform a combined likelihood fit in the full dimen- 
sionality of the parameter space, including variations of important Standard 
Model input parameters. In fact, the feasibility of the approach had already 
been demonstrated in a four-dimensional subspace [8]. This was achieved by 
utilising the Metropolis algorithm in a Markov Chain Monte Carlo (MCMC) 
algorithm [9]. The MCMC approach to parameter sampling has a calculation 
time that scales linearly with the number of dimensions, rendering it very use- 
ful for investigating parameter spaces with dimensionality greater than 3. We 
point the interested reader to refs. [7,9] for details on the algorithm, but the 
important point for the present paper is that the algorithm provides a list (or 
"chain") of points in parameter space. The density of points that it produces 
is proportional to some specified function of the parameters. In reference [7], 
the function chosen was the likelihood, the probability distribution function 
(pdf) of reproducing data d given a point in CMSSM parameter space m: 
C = p(d\m). One ideally wishes to determine the posterior pdf that the model 
point m is correct given the data d: P(m) = p(m\d). The relative posterior be- 
tween two different points m 12 in parameter space can be determined provided 
one assumes a prior pdf for the model parameters p(m): 

Pjmx) _ p(d\m 1 ) p(mi) 
P(m 2 ) p(d\m 2 ) p{m 2 ) 

Thus, the likelihood distribution is proportional to the posterior distribution 
for the case of flat prior distributions, i.e. when p(mi)/p(m 2 ) is a constant 
across the relevant parameter space. For a brief discussion on a Bayesian 
interpretation of the prior probability distributions, see the Appendix. 

Here, we will go beyond the MCMC analysis of Ref. [7] by using a plausi- 
ble natura/ness-favouring prior proportional to 1/c. This will then disfavour 
regions of parameter space with high fine-tuning, arguably a commendable ef- 
fect. We will then compare and contrast the resulting probability distributions 
of observables, parameters and sparticle masses with the flat-prior case. 
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Having obtained a chain of points whose density is proportional to their poste- 
rior assuming some prior, we could in principle just re-weight the probability 
of each point by some new prior and re-plot all of the results. In the present 
case, we would for example use the flat-prior sample, but reweight the prob- 
ability of each point by 1/c, then re-bin any results. The only disadvantage 
to this approach is that if previously unlikely regions of parameter space now 
have high probabilities, they will not have been sampled very often in the 
initial Markov chain run. Thus the statistical fluctuations in this region will 
be large. If the prior doesn't make a lot of difference to the posterior, this 
degradation of the statistics may not matter much. Here though, we will find 
a significant difference in the results due to different priors, implying that a 
simple re-weighting procedure would suffer from poor statistics. We have ac- 
tually re-run the MCMC with the new prior, which circumvents any potential 
statistics problem. 



2 Numerical Methodology 

The numerical analysis follows Ref. [7] identically except for the inclusion 
of a naturalness prior and updates in the codes that calculate observables. 
We briefly include the main features of the analysis here, but we refer the 
interested reader to Ref. [7] for more detail. 

The range of CMSSM parameters considered is sign(/x) = +1, A =-2 TeV to 
2 TeV, m = 60 GeV to 2 TeV, M 1/2 =60 GeV to 2 TeV and tan/3 = 2 to 60. 
Although this range was initially considered on perturbativity and naturalness 
grounds, there is a further pragmatic region for not extending it, namely that 
an enlarged space will result in less efficiency in the MCMC algorithm. The 
MCMC algorithm already uses 4 weeks of CPU time on the lxplus CERN 
cluster, rendering an enlarged run inconvenient. Note that recently, imaginary 
mo has received some attention in the literature [10]. This region predicts 
a non-neutralino dark matter candidate and is not covered by the present 
analysis. 

We use the following measurements in order to constrain the CMSSM: the run- 
ning bottom quark mass in the minimal subtraction scheme m^m),) = 4.2± 
0.2 GeV [3], the top pole mass m t = 172.7±2.9 GeV [11], the strong coupling 
constant in the modified minimal subtraction scheme at Mz a s (Mz) MS = 
0.1187 ± 0.002 [3], the new physics contribution to the anomalous magnetic 
moment of the muon 5(g - 2)^/2 = 19.0 ± 8.4 x 10~ 10 [12,13], the branching 
ratio BR(b -> 57) = 3.52 ± 0.42 x 10" 4 [14] and the relic density of thermal 
dark matter tt DM h 2 = 0. 1126^0^ [15]. 

For each parameter point, S0FTSUSY2.0.4 [16] was used to calculate the 
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sparticle spectrum. The spectrum was linked via the SUSY Les Houches Ac- 
cord [17] to micrOMEGAsl . 3 . 6 [18], which computed predictions for fi^^/i 2 , 
BR{b — > sj) and 5(g — 2)^/2. Empirically derived cuts were then made on 
sparticle and Higgs masses [7] from search data and the requirement of a neu- 
tralino lightest supersymmetric particle. A x 2 value Xi was found for each 
of the six observables mentioned above, and the (un-normalised) posterior 
probability density is calculated as 

P=-J[eM'X 2 m, (5) 

6 i=i 

where c is defined in Eq. 2. The MCMC algorithm was then employed to pro- 
vide a sampling of points in the 7 dimensional parameter space m , My 2 , A , 
tan/?, mb(mb) MS ,m t ,a s (M z )- The density of the resulting parameter space 
points is proportional to P if the chain has converged. 

In order to test the robustness of the numerical results, a convergence test was 
applied. As described in Ref. [7], 9 independent Markov chains of 10 6 potential 
points each were run, all with different random number seeds and at different 
random starting points. An estimation of the convergence of the individual 
chains can be determined by comparing how similar the distributions of inter- 
esting quantities are between the different chains. The technique (invented in 
Ref. [19]) as implemented in Ref. [7], provides an estimated upper bound on 
the potential root variance reduction of each scalar quantity if the chain were 
run for an infinite number of steps. After 10 6 steps, the upper bound is less 
than the required 5% for the naturalness prior sample and for the flat prior 
sample, for every scalar distribution that we are interested in. The chains were 
therefore deemed to have converged. For all of the numerical results presented 
below, the full 9xl0 6 potential-point sample was used. 



3 Posterior Probability Distributions 

In Fig. 1, we show the posterior distribution across various planes in pa- 
rameter space. Unlike ref. [7], we have not shown all possible planes in the 
mSUGRA parameter space, but have instead concentrated on the planes from 
which the effects of the naturalness prior may be deduced. In every plane, the 
posterior probability distribution is shown, marginalised over the other five 
mSUGRA/Standard Model parameters. Marginalisation simply means that 
the unseen dimensions are integrated over. Individual MCMC points were 
placed in 75 by 75 bins in the 2 dimensional space, and for each plot the 
probabilities are normalised by the probability in the highest bin. Thus, one 
can consider the marginalisation process to be an averaging of the posterior in 
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the unseen dimensions. This can be viewed simply as a way of visualising con- 
straints in the 7 dimensional parameter space, or alternatively as a probability 
distribution in the relevant plane, using flat priors in the unseen dimensions. 

On all plots on the left hand side of Fig. 1, the naturalness prior has been used 
in order to calculate the posterior probability distribution. For those on the 
right hand side, a flat prior has been used for the purposes of comparison. The 
results with a flat prior have already appeared in ref. [7] for a previous version 
of SOFTSUSY, and are identical by eye to those shown here. As often pointed 
out in the literature, the CMSSM predicts too much dark matter relic density 
compared to observations unless a specific efficient annihilation mechanism 
acts. 

In Fig. lb, the light Higgs pole region [21,22] is the bright vertical sliver 
at the extreme left hand side of the plot (i.e. low M 1 / 2 and high m ). It 
corresponds to efficient annihilation through an s-channel h°. The stau co- 
annihilation region [23], where fiXi 7T, corresponds roughly to the bright 
region at lowest m , where the lightest stau fi and the Xi are quasi mass 
degenerate. The rest of the probability bulk contains the A — pole region [21] 
where Xi pairs are close in mass to an s— channel A resonance. Some of the 
probability density corresponds to the simultaneous action of several of these 
channels, sometimes including the so-called focus point regime [4] where the 
lightest neutralino contains a significant higgsino component and so Xi pairs 
efficiently annihilate via gauge couplings to pairs of gauge bosons. Slepton 
co-annihilation also contributes in some regions [24]. 

Comparing Figs. la,b, we see that the A°-pole region and stau co-annihilation 
regions have become less probable as a result of the naturalness prior. Also, 
we infer that higher values of My 2 have a penalty originating from the prior 
distribution. Comparing Figs. lc,d, we see the migration of the probability 
distribution away from low mo values toward higher values as a result of the 
naturalness prior. The local maximum in Fig. lc at mo ~ 0.7 TeV corresponds 
to the A — pole annihilation region, which also appears in Fig. Id. Another 
local maximum at m ~ 1.2 TeV corresponds mainly to the increased h°— pole 
region. The third maximum at the lowest values of mo decreases in probability 
as a result of the fine-tuning prior. A comparison of Figs. le,f demonstrates 
again the migration towards lower Mi/ 2 as a result of the naturalness prior. It 
also becomes clear that the /i°-pole region (contained in the horizontal stripe 
at the bottom of both figures) increases in relative importance by using the 
naturalness prior. High tan j3 is somewhat preferred in both cases. 

Although we have discussed the different regions corresponding to various an- 
nihilation mechanisms, the true situation is slightly less clearly defined. In 
order to illustrate this, we use strict definitions of the different annihilation 
regions. We define the h° or A — pole regions by 1 2m x o /m h o/ A o — 1| < 0.1 
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Fig. 1. Comparing constraints on the CMSSM parameter space assuming a natu- 
ralness prior in the left-hand plots and a flat prior in the right-hand plots. Unseen 
CMSSM parameters have been marginalised over. The posterior probability P nor- 
malised to the maximum probability bin of each plot is determined by reference to 
the bar on the right hand side. The contours display the 67% and 95% C.L. contours 
respectively. 
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mechanism 


flat prior natural prior 


h —pole 
,4° -pole 
f — co- annihilation 
rest 


n no£ n c\i 
U.Uzo U.lu 

0.41 0.14 

0.26 0.18 

0.31 0.61 



Table 1 

Relative posterior probabilities for the different co-annihilation regions, as defined 
in the text. The mechanism denoted "rest" indicates the probability of points not 
fitting into one of the clearly defined regions. 

and the stau co-annihilation region is denned by jm^/m^ — 1| < 0.1. Under 
these definitions, the relative probabilities of the annihilation mechanisms are 
shown in Table 1 for each different prior. We see from the table that the nat- 
uralness prior does indeed increase the h°— pole probability at the expense of 
the A°-pole and f co-annihilation regions, but the dominant effect is that most 
of the probability density lies outside the clearly defined regions in the case 
of the naturalness prior. The probability of the h°— pole region is dependent 
upon theoretical errors assumed in the S0FTSUSY2 .0.4 prediction for m^o. The 
LEP2 limits have been softened by 3 GeV in order to take a 3 GeV uncer- 
tainty [25] in the prediction into account in the manner described in Ref. [7]. 
If we instead assume that the S0FTSUSY2 .0.4 prediction of m h o has no the- 
oretical error, the h°-pole region obtains probabilities of 2% and 5% for the 
flat and naturalness prior cases respectively. Investigation of the points in the 
"rest" region shows several competing annihilation mechanisms at work at one 
point: i.e. different combinations of bulk/focus-point/A°//i°-pole annihilation 
mechanisms. For example, the best-fit point in the "rest" region is m = 1419 
GeV, Mi/2 = 180 GeV, A = 550 GeV, tan/3 = 51.8, m b (m b ) = 4.20 GeV, 
m t = 172.4 GeV, a s (Mz) = 0.1188. This point results in a fine-tuning pa- 
rameter of 17.7, BR(b -> S7) = 3.15 x 10" 4 , S(g - 2)J2 = 10.9 x lO" 10 , 
^DMh 2 = 0.1096, indicating a good fit to the data 1 at moderate values of the 
fine-tuning. Investigation of the annihilation mechanism at this point shows 
that annihilation into light quark-anti-quark pairs provides 8% of the anni- 
hilation, bb provides 80% and rf 10%. This is consistent with non-negligible 
contributions from both the h° pole and the Z° pole. The neutralino has a 
non-negligible higgsino component, so the point can also be described as being 
in the focus-point regime. 

Fig. 2 shows probability distributions for the different annihilation mecha- 
nisms. We have fixed the normalisation of each histogram such that its in- 
tegrated probability is 1. The 1/c tendency of each histogram due to the 
naturalness prior is evident at large c. The h°— pole is clearly favoured by the 

1 A good fit is not surprising since we have one more parameter than observables 
used to form the likelihood. 
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Fig. 2. Fine-tuning probability distribution c of various annihilation mechanisms. 
Each individual histogram has been normalised to an integrated probability of 1. 




4000 



Fig. 3. Comparison of mass probability distributions for higgs bosons with ("nat- 
ural prior") and without ("fiat prior") the naturalness prior. Each histogram has 
been normalised to an integrated probability of 1. The plots show distributions for, 
respectively, (a) the lightest CP-even higgs mass and (b) the CP-odd higgs mass. 

naturalness prior, whereas stau co- annihilation and A — poles show higher val- 
ues of the fine-tuning and thus become relatively disfavoured, corroborating 
our interpretation of Fig. 1. A similar investigation of the different compo- 
nents of BR{b — > 57) and (g — 2) M shows much less difference between the 
distributions of the different annihilation regions. Further investigation of the 
MCMC samples finds that 13% of the A pole region is contaminated with 
points that also satisfy the f-co-annihilation constraint. 

Probability distributions for m h o and m^o are shown in Fig. 3. The naturalness 
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Table 2 

Upper 95% C.L. limits on various MSSM particle masses from the fits. Particle 
masses are shown in units of TeV and, except for the lightest CP-even higgs h°, 
have been rounded to the nearest 50 GeV. 

prior skews the h° distribution toward lighter masses. The A distribution 
shape also skews to lower values, except for a small enhanced tail at the 
heaviest masses. Upper 95% C.L. limits on h° and A° masses are displayed in 
Table 2. 



We now turn to the probability distributions for some of the sparticle masses, 
which are shown in Fig. 4. The naturalness prior skews every distribution 
toward somewhat lighter masses except for the right-handed selectron mass, 
shown in Fig. 4d. Along with the other sleptons, this particle becomes heavier 
due to the enhanced probability for the h°— pole region, at large values of mo- 
The squarks in Fig. 4e do not become heavier because their mass is dominated 
by renormalisation group contributions from the gluino mass, which decreases 
with the naturalness prior, as can be seen in Fig. 4c. Enhanced probabilities 
for very low values of m x o and m^± are particularly welcome because they 
enable a future International Linear Collider (ILC) with only 500 GeV centre 
of mass energy to produce them. For a 500 GeV machine, XiXiixtxt) thresh- 
olds can be reached with probabilities of 0.70,0.33 respectively, assuming the 
naturalness prior. For an 800 GeV machine, the probabilities become 0.93,0.58 
respectively. Of course, we should have some idea of the masses of these parti- 
cles from initial data at the Large Hadron Collider (LHC). In order to measure 
sleptons by direct Drell-Yan pair production at the LHC, the sleptons must 
be fairly light: the cross-section a ~ 1 fb -1 for slepton masses of around 400 
GeV would not provide many signal events and the WW fusion rates are neg- 
ligible because of cancellation with Bremsstrahlung processes [26]. The chance 
of directly observing sleptons: m eR < 400 GeV, is 9% (8%) for the flat and 
naturalness priors respectively. From Fig. 4, 95% C.L. upper bounds on the 
relevant sparticle masses are derived as shown in Table 2. Caution must be 
exercised when interpreting upper confidence levels for the scalars: it is possi- 
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Fig. 4. Comparison of mass probability distributions for various sparticles with 
("natural prior") and without ("flat prior") the naturalness prior. Each histogram 
has been normalised to an integrated probability of 1. The plots show distributions 
for, respectively, (a) the lightest neutralino mass, (b) the lightest chargino mass, 
(c) the gluino mass, (d) the right-handed selectron mass, (e) the left-handed squark 
mass and (f) the lightest stop mass. 
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Fig. 5. Comparison of probability distributions of the branching ratio of B s — > [V pT 
with and without a naturalness prior. Each individual histogram has been nor- 
malised to an integrated probability of 1. 

ble that the maximum value taken for mo = 2000 GeV is responsible for the 
upper bound. Presumably, a tail would extend to much higher scalar masses if 
the range of mo were increased. However, the fit itself prefers lighter gaugino 
masses and so they do not have this caveat. From the table, we see that the 
skew towards lighter masses is somewhat moderate. In either case, the bounds 
on the strongly interacting particles imply a good chance of SUSY discovery 
at the LHC [27]. 

In Fig. 5, we compare the probability distribution of BR(B S — > p + p~) as 
predicted by micrOMEGAsl . 3 . 6 for the naturalness and flat priors. The range 
of predicted branching ratios is consistent with the range found by a random 
scan of unconstrained MSSM parameter space in Ref. [29]. The spike at low 
branching ratios is mostly due to the light higgs pole region's contribution 
and is enhanced for the natural prior. The current combined CDF/DO 95% 
C.L. limit [30] 2 BR(B S /x+aO > 3.4 x 10" 7 is displayed by the right-hand 
vertical line, only ruling out a very small amount of parameter space. On the 
other hand, the expected future sensitivity of the Tevatron with 8 fb _1 of 
integrated luminosity is 2 x 10~ 8 [33], shown by the left hand vertical line. 
This would cover 32% of the integrated probability for either the flat or the 
natural priors. 



2 There are newer preliminary CDF/DO bounds [31], for example CDF(DO) have 
non-combined 95% C.L. limits of 2.0(3.0) x 10 -7 respectively. 
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4 Summary 



By employing a simple MCMC, we have investigated the effect of a natural- 
ness prior on CMSSM combined fits. The prior has a significant effect upon 
the fits, indicating that current data are either not sufficiently precise or not 
numerous enough to render posterior probabilities insensitive to the prior. Our 
results bear some similarity to those of ref. [5], which also effectively used a 
naturalness prior, but did not use a likelihood, only LEP 1,2 exclusion limits. 
However, features due to the dominant dark matter annihilation mechanism 
are prominent in our results. The naturalness prior enhances the posterior 
probability of the light CP-even higgs annihilation region of CMSSM param- 
eter space, which comes back inside the 95% C.L. The relative probabilities 
of A — pole and stau co- annihilation regions decrease. Probability distribu- 
tions for sparticle masses are skewed somewhat towards lower values, except 
for the sleptons, which show a slight skew towards heavier masses. The natu- 
ralness prior implies an enhanced probability for the lightest neutralinos and 
charginos to be accessible at a future linear collider. The prospects for evi- 
dence of B s — > at the Tevatron remain good at 32% whichever prior is 
used for 8 fb _1 of integrated luminosity. Rather light CP-even Higgs masses 
are preferred, the 95% upper C.L. limit being only 120 GeV. This is unfortu- 
nately a somewhat difficult regime for h° discovery by the LHC, and implies 
that several years of data-taking will be required [27]. Ref. [28] found that 
fine-tunings of at least 180 are required by the LEP2 bound m h o > 114 GeV, 
tan/5 = 10. Our results indicate fine-tunings in the range of 20-30 are still 
feasible. The results in the present paper include wider variations of more pa- 
rameters (particularly in Aq), but the dominant difference to Ref. [28] is that 
here we allow a 3 GeV theoretical error upon the m^o prediction, which ex- 
tends the fine-tuning range of the fitted points downwards. The random scan 
of Ref. [28] was also very sparse in the parameters and so there is a distinct 
possibility of just not finding the less fine-tuned points. 

The amount of skew of the probability distributions will be modified by chang- 
ing the prior assumption. For example, choosing a prior that is proportional 
to 1/c 2 instead of 1/c (thus effectively counting the naturalness twice rather 
than just once) would produce more skew in the distributions. 

There are many possible future directions for the use of MCMC algorithms 
in fits to supersymmetric models. One obvious direction is the possibility of 
considering different supersymmetry breaking models, for example relaxing 
some of the universality assumed in the CMSSM. Another possible direction 
would be to include electroweak observables in the fit. Yet another would be 
to explore the fi < part of parameter space and work out its probability nor- 
malisation with respect to the \i > part. The MCMC procedure carried out 
here was at the limit of CPU capabilities available to the author, but finding 
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a more efficient algorithm could allow for an enlargement of the parameter 
space considered. In particular, going to high values of m would be useful in 
order to fully investigate the focus point regime [20]. 



Acknowledgements 

We would like to thank G G Ross for providing the idea and motivation for 
this project, to CERN for hosting the computer facilities used and to M Drees 
for an argument about the priors. This work has been partially supported by 
PPARC. 



A Interpretation of the prior 

We are aware that readers that are unfamiliar with Bayesian statistics might 
feel uncomfortable at the appearance of the admittedly subjective priors (see 
for example, a footnote on page 6 of Ref. [32]). A Bayesian viewpoint is that 
we are allowing the data to modify our uncertainty in some parameter and 
that probability is a measurement of our uncertainty. For example, one might 
estimate the chances of throwing a 6 on a straight cubic die to be 1 in 6, but 
knowledge of the position, the speed and the initial configuration of the die 
etc. would allow us to calculate the spin and trajectory of the die and so to 
modify the probability of throwing a 6. 

Solace may be taken from the fact that the more precise data are, the more 
insensitive the posteriors will be to the priors. For instance, the data on the 
electron mass m e is so precise that different priors proportional to, say, ln(m e ), 
e~ me , 1/m.g... will result in almost identical posterior distributions for m e . 
Choosing an extreme prior such as a delta-function for m e at 1 TeV, will 
result in a posterior delta function also at 1 TeV and it appears that many 
feel that this renders the Bayesian approach followed in this paper suspect. 
Usually, the sub-text to this argument is the assumption that the prior is 
some probability distribution set by nature. However, if we interpret the prior 
as a measure of our uncertainty on m e , setting it to be a delta function is 
equivalent to stating that we are absolutely sure of its value and we are not 
going to let data give us any information about it. 
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